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Q_i| Abstract 

, A connection between the algebra of rooted trees used in renormalization theory and 

PsJ ' Runge-Kutta methods is pointed out. Butcher's group and B-series are shown to provide 

a suitable framework for renormalizing a toy model of field theory, following Kreimer's ap- 
, proach. Finally B-series are used to solve a class of non- linear partial differential equations. 

^ ■ 1 Introduction 

' The purpose of this paper is to point out a link between two apparently remote concepts: 

0^ ! renormalization and Runge-Kutta methods. 

0^ ' Renormalization enables us to remove infinities from quantum field theory. Recently, 

Kreimer discovered a Hopf algebra of rooted trees that brings order and beauty in the in- 
tricate combinatorics of renormalization 0|. He established formulas that automate the 
^ il subtraction of infinities to all orders of the perturbation expansion, and proved the effective- 

ly ■ ness of his method for the practical computation of renormalized quantities in joint works 

with Broadhurst |^ and Delbourgo [p|. Moreover, his approach shines new light on the 
problem of overlapping divergences and on the mechanics of the renormalization group 
1^. Furthermore, Connes and Kreimer revealed a deep connection between the algebra of 
rooted trees (ART) and a Hopf algebra of diffeomorphisms . 
I On the other hand, Runge published in 1895 an efficient algorithm to compute the 

solution of ordinary differential equations. For an equation of the type dy/ds = f{y{s)), he 
defines recursively ki = f{yn), ^2 = fiyn + hki/2), j/„+i — yn + hk2- His algorithm was 
improved in 1901 by Kutta, and became known as the Runge-Kutta method. It is now one 
of the most widely used numerical methods. 

In 1972, Butcher published an extraordinary article where he analyzed general Runge- 
Kutta methods on the basis of the ART. He showed that the Runge-Kutta methods form 
a groupie and found explicit expressions for the inverse of a method or the product of two 
methods. He also defined sums over trees that are now called B-series in honour of Butcher. 

Altough the Hopf algebra structure of ART is implicit all along his paper, Butcher did 
not mention itQ. Important developments were made in 1974 by Hairer and Wanner 
Since then, B-series are used routinely in the analysis of Runge-Kutta methods. 

Our main purpose is to show that the results and concepts established by Kreimer fit 
nicely into the Runge-Kutta language, and that the tools developed by Butcher have a range 
of application much wider than the numerical analysis of ordinary differential equations. 

^Hairer and Wanner called it the Butcher group |^]. 

^The opposite of the present Hopf structure of ART was discussed in 1986 by Diir ([0], p. 88-90), together 
with the corresponding Lie algebra (identical, up to a sign, with the one defined in [^). 
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The present expository paper will be reasonably self-contained. After an introduction to 
rooted trees, the genetic relation between ART and differentials is presented. Then Butcher's 
approach to Runge-Kutta methods is sketched. Several B-series are calculated and the 
connection with the Hopf structure of ART is exhibited. The application of Runge-Kutta 
methods to renormalization is exposed using a toy model which is solved non perturbatively. 
Finally, the solution of non-linear partial differential equations is written as a formal B-series. 

2 The rooted trees 

A rooted tree is a graph with a designated vertex called a root such that there is a unique 
path from the root to any other vertex in the tree . Several examples of rooted trees are 
given in the appendix, where the root is the black point and the other vertices are white 
points (the root is at the top of the tree). The length of the unique path from a vertex v to 
the root is called the level number of vertex v. The root has level number 0. For any vertex 
V (except the root), the father of v is the unique vertex v' with an edge common with v and 
a smaller level number. Conversely, z; is a son of v' . A vertex with no sons is a leaf. Rooted 
trees are sometimes called pointed trees or arborescences. 

The tree with one vertex is •, the "tree" with zero vertex is designated by 1. 

2.1 Operations and functions on trees 

An important operation is the merging of trees. If ii, ... ,tk are trees, t = B~^(ti, t2, ■ ■ ■ , tk) 
is defined as the tree obtained by creating a new vertex r and by joining the roots of ti, 
. . . ,tfe to r, which becomes the root of t. This operation is also denoted hy t = [ti,t2, ■ ■ ■ ,tk], 
but we avoid this notation because of the possible confusion with commutators. 

In 0, Connes and Kreimer defined a natural growth operator N on trees: N{t) is the set 
of |t| trees ti, where each ti is a tree with |t| + 1 vertices obtained by attaching an additional 
leaf to a vertex of t. For example iV(l) = • , 

Some trees may appear with multiplicity. 

A number of functions on rooted trees have been defined independently by several authors: 
|t| designates the number of vertices of a tree t (alternative notation is r{t), p{t) and #t). 

Clearly, \B+{ti,t2, . . . ,tk)\ = \ti\ + jtal + • • • + |tfe| + 1. 
The tree factorial t\ is defined recursively as 

• ! = 1 

B+{ti,t2,...,tk)\ = \B+{ti,t2,...,tk)\ti\t2\---tk\. 

(an alternative notation is ^{t)). The notation t\ is taken from Kreimer |^] because t\ general- 
izes the factorial of a number. Besides t\ has also similarities with the product of hooklengthes 
of a Young diagram in the representation theory of the symmetric group [ ]l2[ . A few examples 
may be useful 




2.2 On CM{t) 

CM{t) was defined in |^ as the number of times tree t appears in iV"(l) where n = \t\ is 
the number of vertices oft. In the literature (||l||, [0, p.92, ||l5|, p.l47), CM{t) is written 
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as a{t) and considered as the number of "heap-ordered trees" with shape t, where a heap- 
ordered tree with shape i is a labelhng of each vertex of t (i.e. a bijection between the 
vertices and the set of numbers 0, 1, . . . \t\ — 1) such that the labels decrease along the path 
going from any vertex to the root. This is called a monotonic labelling in [l5[| , p. 147. 




There are {n — 1)1 heap-ordered trees with n vertices. This can be seen by a recursive 
argument. Take a heap-ordered tree t with n vertices, and make n labeled trees by adding 
a new vertex with label n to each vertex of t. Then all the created trees are heap-ordered 
(because the added vertex is a leaf and all other labels are smaller than n). Furthermore, 
all the heap-ordered trees created by this process from the set of heap-ordered trees with n 
vertices are different. Therefore, there are at least nl heap-ordered trees with n + 1 vertices. 
On the other hand, in each heap-ordered tree t with n + 1 vertices, the vertex labeled n is a 
leaf, therefore t can be created from a heap-ordered tree with n vertices by adding this leaf 
with label n. So there are exactly nl heap-ordered trees with n+1 vertices. We shall give a 
non combinatiorial proof of this fact in the sequel. 

Since N{t) is defined by the addition of a leaf to all the vertices of t, a{t) is the number 
of heap-ordered trees with shape t. This number has been calculated in jl^ (see also. [Q, 
p.92): 




where St is the symmetry factor of i, defined in ||, ^ and in where it is denoted by cr{t). 

Note that there is a simple correspondence between the permutations of rt — 1 numbers 
and the heap-ordered trees. Let (pi, . . . ,p„_i) be a permutation of(l,...,rt — 1), then 

• pi is a subroot, labeled pi 

• for i=2 to n — 1 

— if all pj for 1 < j < i are such that pj > pi, then pi is a subroot, labeled pi 

— otherwise, let pj be first number such that pj < pi, in the series pi_i,pi_2, . . . 
then the i-th vertex, labeled pi, is linked to pj by a line 

• when all (pi, . . . ,Pn-i) have been processed, all subroots are linked to a common root, 
labeled 

On the other hand, starting from a heap-ordered tree i, t is arranged so that the set of all 
vertices with a given level number are ordered with labels increasing from right to left. Then 
the permutation is built by gathering the labels through a depth-first search (backtracking) 
of the tree from left to right. For instance, the permutation corresponding to the three 
labeled trees of the above example are (312), (231) and (213). 

Finally, we use the term algebra of rooted trees and not Hopf algebra of rooted trees 
because, thanks to the work of Butcher, the Hopf structure is only one aspect of the ART. 



3 Differentials and rooted trees 

Assume that we want to solve the equation {d/ds)x{s) — F[x{s)], x(sq) — xq, where s is a 
real, x is in M.^ and F is a smooth function from M" to M^, with components f ^{x). This 
is the equation of the flow of a vector field. 
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3.1 Calculation of the n-th derivative 

Let us write the derivatives of the z-th component of x{s) with respect to s: 



^ - i/'N.)i^E|£w»)i^^E|£w»)i/'W.)i 

j j 

jk ^ jk 

A simplified notation is now required. Let 

r = nx{s)] 



so that 



ds ds^ JjJ JjkJ J ^JjJkJ ' 



where summation over indices appearing in lower and upper positions is impUcitly assumed. 
With this notation, we can write the next term as 



A A A 



This relation between differentials and rooted tree was established by Arthur Cayley in 
1857 With this notation, there is a one-to-one relation between a rooted tree with n 
vertices and a term of d"'x{s)/ds^ 

3.2 Elementary differentials 

A little bit more formally, we can follow Butcher (iQ, p. 154.) and call "elementary differ- 
entials" the 5t defined recursively for each rooted tree t by: 

s: - r 

= when t^B+{hM,--- ,tu). (1) 

Using this correspondence between rooted trees and differential expressions, we establish 
the identity: 



d5-. 



N5, = — ^5^^,), 

where N{t) is the natural growth operator of rooted trees defined in Ref. 0. 
So that the solution of the flow equation is 

x{s) — xo + ds' exp[s' N]S, 

J So 

= xo+J2^-^^am{so), (2) 
where \t\ is the number of vertices of t, and a{t) is called CM{t) in [^. 



4 Runge-Kutta methods 



We shall see that sum over trees appear quite naturally with differential equations. So, if 
one is given a function <j) that assigns a value (e.g. a real, a complex, a vector) to each tree 
i, is there a function / such that <j){t) — St- Generally, the answer is no. Consider a function 
(f) such that all components are equal (and denoted also by 0): 

^(•) = 1, 0(J) = a, 0(J) = 6, 

so that for any i, p — 1, /j/^ — a and fjflf'' = b. The first two equations give J2j fj = '^j 
so that the third gives fjflf^ = /j*^ = '^^i ^^"^ 4> cannot be represented as elementary 
differentials (i.e. it cannot be the 5t) of a function J \i h ^ a? . In fact, the number of 
functions reachable as elementary differentials is rather narrow. 

Given such a function (j) over rooted trees, we extend it to a homomorphism of the 
algebra of rooted trees by linearity and (^{tt') — (j){t)(p{t') where the componentwise product 
was used on the right-hand side. If vector flows are not enough to span all possible (p, what 
more general equation can do that? As we shall see now, the answer is the Runge-Kutta 
methods. 

4.1 Butcher's approach to the Runge-Kutta methods 

To solve a flow equation dx{s)/ds — F[x{s)], some efficient numerical algorithms are known 
as Runge-Kutta methods. They are determined by a m x m matrix a and an m-dimensional 
vector &, and at each step a vector Xn is defined as a function of the previous value Xn-i by: 

m 

Xi = Xn-i + h'^aijF{Xj) 

m 

Xn = Xn~i + hy^bjFiXj), 
i=i 

where i range from 1 to m. If the matrix a is such that a^j = if j > i then the method 
is called explicit (because each Xi can be calculated explicitly), otherwise the method is 
implicit. 

In 1963, Butcher showed that the solution of the corresponding equations: 

m 

X^{s) = xo + {s - so)^aijF{Xj{s)) 

m 

x{s) = xo + {s-so)J2bjF{X,{s)), 

3 = 1 

is given by 

X^is) = xo + Y. ^' a{t)tl J2 {t)St (so) 

t I I- j=i 

x{s) = xo-f^^^^^a(t)i!0(t)<5t(so). (3) 

These series over trees are called B-series in the numerical analysis literature, in honour 
of John Butcher (jist, p. 264). The homomorphism (f) is defined recursively as a function of 
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a and b, for i ~ 1, . . . ,m: 

M») = 1 

(j)^{B+{ti---tk)) = ^ Oiji ...ay,0ji(ti)...0j,(ifc) 

m 
1=1 

Comparing Eqs.(^) and it is clear that the Runge-Kutta approximates the solution of 
the original flow equation up to order n if (j){t) ~ l/t\ for all trees with up to n vertices. 

In 1972 jl^, Butcher made further progress. Firstly he showed that Runge-Kutta meth- 
ods are "dense" in the space of rooted tree homomorphisms. More precisely, he showed that 
given any finite set of trees Tq and any function 9 from To to M, then there is a Runge-Kutta 
method (i.e. a matrix a and a vector b) such that the corresponding (f> agrees with 9 on Tq 
(see also [Q p. 167). 

4.2 Further developments 

Furthermore, Butcher proved that the combinatorics he used to study Runge-Kutta methods 
in 1963 was hiding an algebra. If (a, 6) and {a' ,b') are two Runge-Kutta methods, with 
the corresponding homomorphisms </> and 0', then the product homomorphism is defined (in 
Hopf algebra terms) by 

{4)i<(t>'){t) = TO[(<?!>0 0')A(t)]. 

Butcher proved that the derived from Runge-Kutta methods form a group. Again, 
this is nicely interpreted within the Hopf structure of the ART. For instance, the inverse of 
the element is simply defined by (f>^^{t) — (j)[S{t)], where S is the antipode. This concept 
of inverse is quite important in practice since it is involved in the concept of self-adjoint 
Runge-Kutta methods, which have long-term stability in time-reversal symmetric problems 
([|l|, p.219). The adjoint is defined within our approach by (/>*(i) = (-l)l*l(/)[S'(i)]. 

On the other hand. Butcher found an explicit expression for all the Hopf operations of 
the ART. Given the method {a,b) for </>, he expressed the method {a' ,b') for (p o S {S is the 
antipode) in (simple) terms of (a, 6). Moreover, (Q, p. 312 et sq.), if (a, 6) and (a', 6') are 
two Runge-Kutta methods (with dimensions m and m', respectively), corresponding to </> 
and (/)', the method (a",&") corresponding to the convolution product (</> ★</)') is 



< 


— Gij 


if 


1 < i < m and 


1 < J < TO, 






= <j 


if 


TO+l<i<mH 


- to' 


and TO + 


1<J< 


< 




if 


TO+l<i<TO + 


to' 


and I < j 


< TO, 


n 


= 


if 


1 < i < m and 


to ^ 


- 1 < j < m - 


h TO , 


h'l 


= b. 


if 


I < i < m, 








h'l 


= K 


if 


m + l<i<m + 


m'. 







In 1974, Hairer and Wanner ([T^], p. 267) built upon the work of Butcher and proved the 
following important result: if we denote 

B{cf>,F) = l + Y^^-l^^am<p{t)St{so) (4) 

then 

5(0', 5(0, F)) = i?(0*0',f^), (5) 

where B{(j)' , B{(t>, F)) is the same as Eq.(|4|), with (t){t) replaced by (t)'{t) and 6t replaced by 
S'{t) (i.e. 5'{t) is calculate as 6t, but with the function _B(0, F)(s) instead of the function 
Fix{s))). 

In other words, the group of homomorphisms acts on the right on the functions F. 
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5 The continuous limit 

In his seminal article [0, Butcher did not restrict his treatment to finite sets of indices. It 
is possible to consider the continuous limit of Runge-Kutta methods. A possible form of it 
is an integral equation, which we write artitrarily between and 1: 

Xu{s) = xo + {s-so) / dva{u,v)F{Xy{s)) 

Jo 

x{s) = xq + {s~so) / b{u)duF{Xu{s)), 
Jo 

the solution of which are 

Xu{s) = a;o + ^ ^'^ J.f a{t)t\ I dva{u,v)(l)y{t)6t{so) 
f r\- Jo 

{s — So)'*' 

x{s) = ^0 + ^ a{t)t\(j){t)St{so). 

The homomorphism is defined recursively as a function of a and b: 
(bui») = 1 

• ■ -ifc)) = / duia{u,ui)(j)ui{ti)- ■ ■ I duka{u,Uk)(l)uk{ik) 
Jo Jo 

(l){t) = f dub{u)(j)u{t). 

Jo 

Continuous RK-methods do not seem to have been much used, except for an example in 
Butcher's book (Q p.325). 

5.1 Butcher's example 

It will be useful in the following to have the results of a modified version of Butcher's example. 
So, we consider: 



Xu{s) - xo + (s-so)/ F[X,{.s)]dv (6) 

Jo 

x{s) = xq + (s - So) / F[Xu{s)]du, 



which corresponds to a(u,v) = l[o.u]{v), b{u) = 1. This Runge-Kutta method will be used 
again in the sequel, and will be referred to as the "simple integral method" . 
If we take the derivative of Eq. @ with respect to u we obtain 

-^X„(s) = (s-so)F[X„(s)], 
du 

so Xu{s) = y{so + (s — so)u), where y{s) is the solution 

y(s) ^xo+ f F[yis')]ds'. 

J So 

Moreover 

x{s) = a::o + (s - so) / F[Xu{s)]du 
Jo 

= xo + {s- So) / -F[?;(so + (s - so)u)]du 
Jo 

= Xo+ f F[y(s')]ds' = y(s). 

J So 



The corresponding homomorphism (j){t) is defined by 



<t>{t) = / du4)u{t). 
Jo 

Using the facts that \B+{ti ■ ■■tk)\ = {\ti\ H + \tk\ + 1) and B+{ti ■ • -tfe)! = {\ti\ H + 

\tk \ + 1)^1 ! • • ■ ifc! it is proved that the solutions of these equations are 

9u[t) = 

If we introduce (j){t) = 1/tl into Eq.(||) we obtain Eq.(||). So we confirm that the sohition 
of the equation 

x{s) =xq+ f F[x{s')]ds' 

J So 

is 

x{s) = XQ+^^^—0^a{t)St{so). 
5.2 First applications 

The above example can already bring some interesting applications. But we must start by 
giving a way to calculate (5t(so) in a simple case. 

5.2.1 Calculation of 6t{so) 

To obtain specific results, we must choose a particular function F. The simplest choice is to 
take a vector function F, with components p{x) — f{J2j ^i/-^)) where N is the dimension 
of the vector space and / has the series expansion 

^ /(")(0)g" 
J[s) = V j ■ 

From the definition of St in Eq.(|l]), one can show recursively that, for i = 1, . . . , TV, 51(0) 
is independent of i (and will be denoted 6t) and 

S, = /(O) 

St - f<^^\0)St,5t2---K when t = ts, • • • , ife)- (7) 

In ref. j^, Kreimer defined a similar quantity, that he called Bt- Here St and Bt will be used 
as synonymous. 

The simplest case is f{s) = exps and sq = 0, where /"-"'(O) = 1 and St = I for all trees t. 

5.2.2 Weighted sum of rooted trees 



If we take / = exp, sq = and xq = in Butcher's example (see section 5.1), we have to 
solve the equation 



!(s) = / exp[a;(s')]ds' 

^0 



which can be differentiated to give x'{s) = exp(a;(s)) with a::(0) — 0. This has the sohition 

x{s) = -log(l - s) = ^ 



n 

n=l 



On the other hand, the corresponding homomorphism is (j){t) = l/t\ and the B-series for this 
problem is 



s 

t 

Comparing the last two results, we find 



\t\=n 

in other words, the number of heap-ordered trees with n vertices is (n — 1)!. 
5.2.3 Derivative of inverse functions 

We can try to extend the last example to an arbitrary function f{x). The equation to solve 
becomes 



x{s) = / f[x{s')]ds', 
Jo 

or x'{s) = f{x{s)) with a;(0) = 0. Let 

'^'^ = Jo wr 

which gives us s = S(x), or x{s) — S^^{s), where is the inverse function of S. If f — exp, 
S{x) = 1 — exp(— cc) and we confirm that x{s) — — log(l — s). 

We can use this result to calculate the derivatives of a function x(s), given as the inverse 
of a function S{x). To do this, we define f{x) = 1/S'{x) and, using Eq.(^), we obtain 

x(")(0) - 5]a(t)<5t, (8) 

\t\=n 



where 5t is calculated from /(s) using Eq.(|^) in section 5.2.1 



This method can also be calculated to find the function / satisfying given values for 



a„ = ^ a{t)dt 

\t\=n 



where St is calculated from /. For instance, if we want 



a(t)St = n!, 

|t|=" 



we must take /(s) = 1 + 



5.2.4 Other sums over trees 

We give now further examples of sums over trees, that will be used in the sequel. For instance, 
assume that we need to compute 



\t\=n 
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This term comes in the Butcher series with (f>{t) — l/{t\)'^. Since this (p{t) is the square of the 
previous one, the corresponding Runge-Kutta me thod can be reahzed as the tensor product 
of two "simple integral methods" (see section In other words 



a{u,u',v,v') = a{u,v)a{u' ,v') — l[Qu]{v)l[Q.u']{v') b{u,u') — b{u)b{u') ^ 1. 
and the Runge-Kutta method is now 

Xuu'{s) = xo + is-sn) dv dv' f[Xyy,{s)] 

Jo Jo 

x{s) = xo + (s - So) / du duf[Xuu'{s)]. 
Jo Jo 

The corresponding homomorphism (p{t) is given by 

(t>uu' ( • ) = 1 

(t>uu'{B^{ti- ■ -tk)) = dui du'^(j}u^u'^{ti)- ■ ■ j duk / du'^(l)u^u>^{tk) 

Jo Jo Jo Jo 

4>{t) = I du I du'4>uu'{t)- 



The solutions of these equations are 

(pun' (t) - — 

so that, from Eq.(^) 

Xuu'{s) = |t|! 1\ 

The conclusion is that X„„/(s) is in fact a function of uu' and not of u and u'. More 
precisely, we know from the general formula Eq.(H) that the B-series for the solution of 

x{s) = xo + {s- So) Jo du /g du' f[Xuu'{s)] is 

/ N - So)'*' a(t) . , . 
x{s) = xo+Y.-^ ^-^((^o), 

so that X„t('(s) = xi^SQ + (s — so)uu'). If we use the successive changes of variables w ~ uu', 
v' = So + (s — so)w and w = sq + (s — so)u we find 

x{s) = + (s — So) I du I du' f[x{so + (s — so)uu')] 
Jo Jo 

= xo + {s ~ So) — / du;/[a::(so + (s - So)w)] 







u 



rill rso + {s-sa)u 

X0+ I — dv'f[x{v')] 







u 



= X0+ r f" dv'f[x{v')]. 

J so « - So J so 

With the initial values xq ~ sq ~ this gives us 



xis) =1^1 dv'f[x{v')], (9) 
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or sx" + x' = f{x) with a;(0) = and a;'(0) = /(O). If we take again f{x) — exp(x) wc find 
sx" + x' = exp(a;) with a;(0) = and x'{0) ~ 1, so that 

oo „ 

x(s) = -21og(l~s/2) = ^ 
Comparing this with the B-series 



n=l 



X{S 



^ a{t) 
^\t\\ t\ 

we obtain 



\t\=n 



which is the result found by Kreinier in using combinatorial arguments. 

As a final example, we can consider the Runge-Kutta method a{u,v) = 1, b(u) = 1 
which gives (j){t) = 1 for all trees t. The equation for x{s) is now a fixed point problem 
x{s) = sexp(a;(s)), whose well-known solution is 



s 

n" 

n\ 



so that 

|t|=rj 

These examples show that B-series can be used as generating series for sums over trees. 



5.3 The antipode 

The Hopf algebra structure of the ART entails an antipode S. If <^{£) is an homomorphism, 
the action of the antipode on can be written as S{(t)){t) = (j){S{t)). If the Runge-Kutta 
method for </) is B, then the Runge-Kutta method for (p^ — S{(j)) is = Au — B, 
B^ = —B. It is useful to see it working on simple cases: 

^uil) = ^f(<^f(«))-A„(0.(.))-S(0,(.)) = 04 J) -</.(•) 
^'il) = -i3(0f(I)) = -0(J) + S(l)0(.) = -0(J) + 0(. )</>(.) 

5.4 The convolution 

The convolution of 4> and (/>' is defined as 4'"{t) = {4>* 4>'){i) = rn[{(l) (g) 0') A(t)]. 

Let Au, B and A'^, B' be the Runge-Kutta methods of, respectively, and 4)'{t). To 
be specific, we consider that u varies from to 1. Then the Runge-Kutta method for 0" is 
A", B" , where u varies from to 2 and 

AuiXv) if 0<u<l and 0<u<l 
if < It < 1 and 1 < w < 2 
B{Xy) if 1 < u < 2 and < w < 1 
A'u_i{X^^i) if l<u<2 and l<w<2 
B{X^) if < u < 1 
B'{X^^i) if 1 < w < 2. 



K{X.) = 

K{X.) = 

K{X.) = 

B"{X,) = 

B"(X,) = 
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Again, we show the formula in action: 



0"(.) = + + 

n]) = B{A^{i))+B'mi) + K_,m = b{m*)) + b{i)b'{\)+b'{^'u{*)) 
= 0(J) + 0(.)0'(.) + 0'(I). 

6 Runge-Kutta methods for renormalization 

In this section, we shall follow closely Kreimer's paper and define, for each operation on 
homomorphisms, a corresponding transformation of the Runge-Kutta methods. Instead of 
attempting a general theory, we consider a specific example in detail. 

6.1 Runge-Kutta method for bare quantities 

We consider that a given bare physical quantity can be calculated as a sum over trees, and 
that the corresponding Runge-Kutta method has been found as a pair of linear operators 
Au and B. The usual combinatorial proof show that the solution of the equations (we take 
so = 0) 

X^{s) = xo + sA^[f{X,{s))] 
x{s) = xo + sB[!{X^{s))l 



IS 



Xu{s) = XQ+^—a{t)t\Au[(t>v{t)]5t 
t I'^l- 

x{s) = xo + }_^^a{t)tmt)5t, 



where, as usually. 



0n(«) = 1 

(j3u{B+{ti---tk)) - Au[(Pm{ti)] . . . Au[(l}u,{tk)] 

<P{i) = B[<P^{t)]. 

Here x{s) is the sum giving the bare quantity of interest. In the examples developed by 
Broadhurst and Kreimer ||^, the quantity of interest is 



where Bt is obtained recursively from given _B„ by 
B, = Bi 

Bt = B^t^St,St2 ■ ■ ■ St, when t = B+ {ti,t2, ■ ■ ■ (10) 

In the renormalization problems considered by Broadhurst and Kreimer, the _B„ are 
defined from a function L{d) regular (and equal to 1) at the origin, by 

On — ■ 

ne 
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A pair of operators giving — Bt can be defined as 



e Jq dv 



The quantity of interest x{s) is then obtained by tensoring Au with the "simple integral 
method" to obtain (j){t) = Bt/t\. 

The only thing that we need in the following is the action of A„ on a monomial v"^^ 

A„(t;"-i) ^ - f L(e-^w)v"-i = - f v''-^dvL{ne) = B„m". (11) 
e Jo "u e Jo 

6.2 Sr, the "renormalized antipode" 

In ref. ^ , Kreimer defines recursively a renormalized antipode^ depending on a renormaliza- 
tion scheme R. We take as an example the toy model used by Kreimer, where R[4>\ — {(j)) is 
the projection of on the pole part of the Laurent series in e inside the bracket. 

Following the results of section |3.3| , the Runge-Kutta method for Sf/icj)) can be obtained 
from the Runge-Kutta method of by A^{X) = Au{X) - {Ai{X)), B^(X) = -{Ai{X)). 
Working out the first examples using Eg. (pi]), we find, 

^^•) - 1 

^^(.) = (Ai(l)) = (Bi) 



(J) - AM»))-{Ai{M»)))=Au{l)-{Ai{l))^B,u-{B,) 
(J) = -(Ai(</)f(J)))=-(i?2i?i) + ((i?i)i3i). 



6.3 Renormalized quantities 

Finally, the renormalized quantities x^{s) are obtained from the convolution of Sfj{(l>) with 
(j). To obtain the corresponding Runge-Kutta method, we use the results of section ^.4| . 
However, the domain where 1 < u < 2 is not used, and the Runge-Kutta method for the 
renormalized quantity is A^{X) = A^{X) ~ {Ai{X)), B^{X) = Ai{X) - {Ai{X)). It may 
seem surprising that such a simple equation encodes the full combinatorial complexity of 
renormalization. It is not even necessary to work examples out, because A^{X) — Af^{X) 
so that 0^(i) = (/>f (i), and the only difference comes from the action of B^. 

For a real calculation of a;^(s), we do not need A^ and B^ which give us 0(t) = r(t), but 
the tensor product of this method with the "simple integral method" to obtain 0(t) = tW{t). 
In detail, the equation for the renormalized quantity x^{s) is 

X^„,(s) = - r dv r dv'L{ed,v)e^^^'^'^ - {- f dv C dv' L{td.,v)e^^^'^'^) (12) 
^ Jq Jo £ Jo Jo 

x^is) = - f dv f dv'L{ed.,v)e^^^'^''^ - {- [ dvf dv'L(ea„w)e^-'(^)). (13) 
^ Jo Jo ^ Jo Jo 



For a general renormalization scheme R, one replaces {Au{X)) by R[Au{X)]. Finally, 
Chen's lemma for renormalization schemes |^ is obtained from Hairer and Wanner's theorem 
Eq.(D. 



7 Renormalization of Kreimer's toy model 

In this section, we use Runge-Kutta methods to renormalized explicitly Kreimer's toy model 
for even functions L{e). In remarkable properties of the renormalized sum of diagrams 
with "Connes-Moscovici weights" were noticed. 

^In Hopf algebra terms SR{(j>){t) = -R[<t>{t) + m[{SR (g) /d)((?!) » 0)P2A(f)]. 
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7.1 Equation for the renormalized quantity 

The role of the sum over u' m Eg. (p^) is to add a factor l/t\, as in section ^^2!4 . Therefore, 
the same reasoning can be used to show that X^^, (s) is in fact a function of su' and we 
write X^^,{s) = X^{su'), which defines the function X^{s). The equation for X^{s) can 
be found from Eq.Q and the relation X^(s) = ^^^(l) as 



X„ 



\s) = - f dv f ds'L(ea„w)e^"(''') - (- / dv f ds' L{ed^v)e^-^'">) (14) 
^ Jo Jo ^ Jo Jo 

\s) = - f\v r ds'L{edyv)e^-^'"> ^ {- f\v r ds'L{ed,v)e^-^''^). (15) 
^ Jo Jo ^ Jo Jo 

To solve this equation, we expand X^{s) in a power series over u: 

oo 
n=0 

A standard identity gives us 

oo 

exp(X^(s)) ^A„(a)M", where 

A„(a) = ^ with \a\ = ai + 2a2 -\ + na„. 

^-^ ai! • • ■ a„! 

\a\—n 

A„(a) depends on s through its arguments ai(s). The sets of ai for a given n can be obtained 
from the partitions of n: . . . , /i„), where /ii > • • ■ > /^n by a„ = /x„, ai — [li — fii+i for 
i < n. 

The first A„(a) are 

Ao(a) = 1 Ai(a) = ai A2(a) = 02 + y hia) = 0,3 + aia2 + ■ 

7.2 Solution of the equation 

Introducing the series expansions for X^{s) and exp(X^(s)) into Eq.(|l|) we obtain 

00 00 ^5 00 

^a„(s)?/" = / e-"^'">K{a)ds'u^+'-{Y,Bn+i e^"^'") Xn{a)ds') 



n=0 

or 



ao(s) 



a„(s) = B„ / e''"("')A„_i(a)ds' for n > 0. (16) 





To solve this equation, we need to go back to the equation for the bare quantity 

X°(s) = - r dv r ds'Lied,v)e^°^''\ (17) 
f Jo Jo 

Again ^^(s) is a function of su, we define X'^{s) = (1) which satisfies 

^ Jo u Jo 
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The solution of this equation is given by the B-series 

X'{s) ^ 5]a„." with a„=^^i|^. (18) 

On the other hand, we can also expand e'^"*^^' using the functions A„(a). Identifying both 
sides of Eq.([l8|), we obtain the relation 

a„ = — A„-i(a). (19) 
n 

With this identity, we can now prove that, for the renormalized quantities, 

a„(s) = {g{s))"-an, where g{s) = / exp{ao{s))ds' . (20) 

Jo 

Since Ao(a) — 1 and cin ~ Bi, this equation is true for n — 1, from Eq. (p^ . If Eq.(pO|) is 
true up to n — 1, then A„_i(a) = (.g(s))"^^A„_i(a) and the derivative of Eq.(p^) gives us 

a'^is) = i?„e'^°(^)A„_i(a) = B„5'(s)(.g(s))"-iA„_i(a) = n(g(s))"-ia„, 

by Eq.(p!9|). Integrating this equation with the condition a„(s) = gives Eq.(^ at level n. 

By this we have proved that the flow for the renormalized quantity is a reparametrization 
of the flow for the bare quantity: X^{s) = ao{s) + X° {ug{s)) and X-'^(s) = ao{s) + X°{g{s)). 

To determine ao(s) we proceed step by step. In Eq.(|l8|) we expand L{edyv) over e. The 
first term is just 1, and we obtain Eq.(^) with the solution x{s) = —2 log(l — s/ (2e)). For the 
renormalized quantity, the most singular term becomes X^{g{s)) — — 21og(l — g{s)/{2e)). 
Since X^{s) is regular, this singular term must be compensated by a corresponding term in 
ao(s). By equating the most singular terms we obtain ao(s) = — 21og(l — g{s)/{2e)). We 
know from Eq.(po|) that ao(s) ~ log(g'(s)), and we obtain the most singular terms as the 
solution of g'{s) = 1/(1 — g(s)/(2e))^, which is: 

9{s) = ^ 

ao(s) = -21og(l + ^). 

By expanding ao(s) as a series in s, we obtain the most singular term observed in |^ and 
proved in j^. One notices that the singularity of the non-pertubative term ao(s) is logarith- 
mic, and much smoother than the singularities coming from the expansion over s (i.e. the 
perturbative expression). 

7.3 Differential equation for the finite part 

In general, one should proceed now with the next singular term. To obtain it we denote 
Y{s) — X^{g{s)), this change of variable gives the equation for Y{s): 

Yis) = ir^rd.,'(.)x(.+.4^a.)e-(^). 

e Jo 9iu) Jo .9 («) 

Now we can write Y{s) = X^{s)~ao{s), and notice that the term — ao(s) on the left-hand 
side is compensated by a term on the right-hand side where L — 1 and exp(X^(s)) = 1. We 
obtain the equation for X^{s): 



1 .,.a.„4a„„l.i;)%-"... 



The nice aspect of the previous equation is that it seems to have a limit as e goes to zero. 
In fact, it has a limit when L is even, as we shall show now. 
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Writing X{s) — lim^^Q X^'{s), and taking the limit e — > in the previous equation, we 
obtain 



Xis) 

or, in differential form: 



2 / ^ 



dv 



-As'X'is))' = 



s2^^ 2 ds 



If X{s) and L{S) are expanded as 



X{s) = ^6„s" and L((5) = 1 + ^ i„5", sothat 



ra=l 



= 1 



2 ds 



(21) 



we obtain the following relation for the term in s: bis — (bi + Li/2)s. If Li is not zero, we 
obtain a contradition and must proceed with the withdrawal of divergences. For simplicity, 
we shall assume that Li = 0. Then bi becomes a free parameter of X(s). All terms bn with 
n > 1 can now be determined from bi and L„ (n > 1). All terms are regular. 

In [Q, the function L(S) was taken even. Then Li — 0, and their results correspond to 
bi — 0. Broadhurst and Kreimer have also used a function L(e,S). The present treatment 
can be applied to this more general situation, with the only change that 

Lri = n\ lim lim — ^L(e, S). 

Clearly, Eq.(pl|) is much faster to solve than computing the sum over trees. For instance, 
the expansion could be calculated up to 20 loops (i.e. 620) within a few seconds with a 
computer. 



7.4 Alternative point of view 

There is an alternative way to solve Eq.(p7|) for the bare quantity. We define a function /(s) 
from L{5) by 

n=0 

A relation between /(s) and L{S) can also be established through the Mellin transforms of 
/ and L as M{f){z) = M{L){t - ez)T{z). 

With /(s) we can write the equation for the bare quantity as 

X\s) = 1 fdvfiX^iv)). (22) 

Alternatively, one can go from f to L and consider the results of the toy model as a 
method to renormalize equations of the type (p2[). 



8 n-dimensional problems 

For applications to classical field theory, we need to develop Runge-Kutta methods for the 
n-dimcnsional analogue of the flow equation: non-linear partial differential equations. The 
purpose of the present section is to indicate how B-series can be used for this case^. The 
method apply to equations of the form Lipir) = F['0(r)], where L is a differential operator 
(e.g. the nonlinear Schrodinger equation A?/; = 4''^)- 

''Kreimer was independently aware of the possibility to use B-series for non-linear partial differential equations. 
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8.1 Formulation 



We need two starting elements: a function V'o(r) which is the solution of i'i/;o(r) = 0, and a 
Green function G(r, r'), that is a solution of the equation LrG{r, r') = S{r — r'), with given 
boundary conditions. The function V'o(r) will play the role of an initial value, and the Green 
function will decide in which "direction" you move from the initial value. It will also state, 
in some sense, the boundary conditions of the solution "(/"(r). 

Using these two functions, the differential equation L'ip{v) = F['i/;(r)] is transformed into 
■j/;(r) = ■i/'o(r) + / (ir'G'(r, r')F[-i/;(r')]. The action of L enables us to go from the second to 
the first equation. 

The combinatorics is the same as for the standard Runge-Kutta method, and the result 

is 



V'(r) = Vo(r) + }_4;^ / dr'G(r,r')(/..'(t), (23) 
where (j)r{t) is defined recursively by 

MB+iti---tk)) = F'^''\Mr)]J dnG{r,ri)(br,{tiy-- J drkG{r,rk)(br,{tk). (24) 
If -0 is a vector field, the solution is the same, and equations (|2j) get indices: 

0;(«) ^ fiMr)] 

ct>UB+{h ■ ■ ■ tk)) = [^o(r)] / dnG'^ (r, r,)<pi\ (h)- ■ ■ ( dr^G^ (r, r,)^i (tk), 



where (r, r') is a component of the matrix Green function. 

In the previous sections, the series (||) was written as a function of (j)(t) (describing the 
effect of the Runge-Kutta method (a, 6)) and St (describing the effect of the function F[a;]). 
In the present case, this separation is no longer possible, and 4>{t) combines both pieces of 
information. 

8.2 Examples 

In this section, equation (|2^) is applied to the one-dimensional problem and to the Schrodinger 
equation. 



8.2.1 The one-dimensional case 

It is instructive to observe how the one-dimensional case is obtained from Eq.([23|). The 
differential operator is L = d/ds, so the initial function tpois) must satisfy d/ds^pois) — 0: 
V'o(s) is a constant that we write xq. For the Green function G(s, s'), we have the equation 
LG{s,s') — 5{s — s'), so G{s,s') — 9{s ~ s') + G(s'), where 9{s) is the step function and 
G(s') a function of s' . To determine G(s'), we note that, in the "simple integral method", 
there is an integral from sq to s. From the Green function G(s, s') — 0{s — s') — 6{so — s'), 
we obtain 

/OO PS 
G{s,s')fis')ds' = / f{s')ds' 
-OO Sq 

which is the required expression. 

Now, the role of ipo and the Green function is clear for the one-dimensional case: ipo gives 
the initial value xq and G specifies (among other things) the starting point sq- To complete 
the derivation of the one-dimensional case, we note that ipo{s) = xq does not depend on s, so 
the terms F^^) [V'o(s)] = F('=)[ Xq] are independent of s and can be grouped together to build 
St as in (||). On the other hand, the integration over Green functions build up (s — sq)'*' /t\ 
and we obtain Eq.(||). 
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8.2.2 The Schrodinger equation I 

If we write the Schrodinger equation as {E + A)V'(r) — V(r)^{r), we can apply Eq.(|25 
with F[ip] = V{r)ip. We take for 0o(r) a solution of {E + A)V'o(r) = and for G(r, r') the 
scattering Green function (e.g. G(r — r') = —e*^!'"^'' l/(47r|r, r'|) in three dimensions). 

The calculation of (f>{t) is straightforward because, in a such a linear problem, F^''^ = 
for /c > 1. Hence, the only rooted trees that survive are those with one branch. For these 
trees a{t) — 1 and t\ — \t\l and we obtain 



^(r) = ^o(r)+ J driG(r,ri)y(ri)V;o(ri) + J dridr2G(r, ri)y(ri)G(ri , r2)y(r2)^o(r2) + ■ • 
where we recognize the Born expansion of the Lippmann-Schwinger equation. 
8.2.3 The Schrodinger equation II 

We can also treat the Schrodinger equation in an alternative way as the system of equations: 

{E + A)ij{r) V{p)i^{r) 
dpi 

This is a matrix differential equation. We give index to the first line, and index i (running 
from 1 to the dimension of space) to the other lines, called the space lines. The purpose of 
the space lines is just to ensure that p = r. This is a standard trick to take the r dependence 
of V into account in the expansion (see e.g. ||l5| p. 143). As initial value we take "00 (r) and 
Po = 0, the matrix Green function is diagonal and it is equal to the scattering wave function 
for line and to 6{ri — r-) — S{—r[) for line i. 

For (/!),■(•), the zero-th component is y(0)V'o(r) and the space components are 1, for all 
the other trees, the space components are and the zero-th component of the simplest tree 
is 

Ml) = V{Of j dv'G{v,v')i^^{r') + Y,rAV{Q)M^) 

i 

M^) = V{Of I driG{r,ri) + V{0) I dr'G{r,r')Y,rAV{0)Mr') 
MJ\) = 2V{Q)Y,nd,V{Q) j dv'G{v,v')Mr')+Y.nr,d,d,V{Q). 

i ij 

The expressions become more and more complex, but their derivation is made systematic by 
the recurrence relation. 



9 Conclusion 

Butcher's approach to Runge-Kutta methods was applied to some simple renormalization 
problems. Since Cayley, it is clear that the ART is ideally suited to treat differentials. This 
was confirmed here by presenting a B-series solution of a class of non-linear partial differential 
equations. 

The recursive nature of B-series make them computationally efficient: 4>u{t) can be ob- 
tained by a simple operation from the (j)u{t') of smaller order t' . This is why B-series can be 
automated and implemented in a computer. 

Butcher's approach has still much to offer. In the numerical analysis literature, B-series 
have been generalized to treat flow equations on Lie groups. The main change Q is to 
replace the algebra of rooted trees by the algebra of planar trees (also called ordered trees 
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|l9|). The elementary differentials get then a "quantized calculus" flavor, especially in the 
definition given Munthe-Kaas in terms of commutators with the vector field F = pdi (see 
also Ginocchio). Using this generalized ART, extended work has been carried out recently 
for the numerical solution of differential equation s on Lie groups (see Ref. [|l8|, ^ and the 
web site jittp : / / www . math . ntnu . no/num/ synode ) . 



B-series have been generalized in other directions, e.g. stochastic differential equations 
^ and differential equations of the type dy/ds = f{y,z), g{y,z) = 0, which are called 



differential algebraic equations |22 



It is our hope that Butcher's approach can be applied to quantum field theory. 
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11 Appendix 

For further reference, the action of the coproduct and the antipode on the first few trees are 
given here. 

11.1 Coproduct 

Al = 1 (g) 1 
A« = •®l + l(g)« 
aJ = J®l + l(8)J + «®« 

A^ = ^ ® 1 + 1 (g) ^ + J ® • + • (g) J 
^<A) " ^® 1 + 1 • • ® • +2« g) J 

a| = I g) 1 + 1 g) I + • g) ^ + ^ g) • + J (g J 

^ ^^gi 1 + 1 g)^^+ 2« g) ^ +^g) • + • • g) J 
A^^ = ^\g) 1 + 1 g)^\+ • g) ^ + • g)^+ J g) J + • • g) J + • J g) • 
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11.2 Antipode 



I 



51 



^+2. J 



• • • • 



S 



• • • • 



• • • • 
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